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ROBUST  MULTIOBJECTIVE  CONTROLLERS  FOR 
SMART  STRUCTURES 


Vittal  S.  Rao,'  Robert  Butler,’  Wei  Zhao’ 


ABSTRACT  Flexible  smart  structures  are  large  mechanical  structures  with 
applications  in  which  specific  performance  characteristics  are  desired  in  the  presence  of 
parameter  variations  and  disturbances.  These  structures  tend  to  have  severe  controller 
effort  restrictions  and  tightly  spaced  lightly  damped  modes.  When  designing  controllers 
for  smart  structures,  H2  controllers  are  well  suited  for  control  effort  restrictions  and 
closed  loop  performance  specifications  while  controllers  are  better  suited  for 
modeling  uncertainties.  This  paper  examines  the  application  of  H2/Ho<,  controller  design 
methodologies  to  smart  structures. 

A  unique  feature  of  mechanical  distributed  systems  is  that  state  space  systems  can  be 
determined  from  finite  element  models  (FEM)  in  which  the  states  have  physical 
significance.  For  certain  systems,  these  states  can  be  directly  measured  by  using  a 
distributed  PVDF  film  appropriately  shaped  and  applied  to  the  structure.  This  full  state 
feedback  system  allows  for  the  implementation  of  H2/H0.,  full  state  feedback  algorithms. 
The  development  of  a  control  system  implementing  this  type  of  algorithm  is  described 
for  a  simple  cantilever  beam. 

In  addition,  a  H2/H^  algorithm  developed  by  Bernstein  and  Haddad  is  investigated 
for  the  cantilever  beam.  This  method  does  not  require  state  measurement,  but  the 
controller  design  algorithm  is  slightly  more  complicated.  This  algorithm  requires  the 
solution  of  one  algebraic  Riccati  equation  and  two  coupled  algebraic  Riccati  equations 
which  are  solved  iteratively. 

INTRODUCTION 

Control  of  smart  structures  presents  a  number  of  difficult  problems  which  arise  due 
to  structural  parameter  variations,  limited  actuation  force,  nonavailability  of  states  for 
feedback  and  restrictive  control  bandwidth.  Several  robust  control  system  design 
techniques  have  been  developed  and  utilized  in  the  design  of  controllers  for  smart 
structures  (Lashlee  et  al.,  1994;  Damle  et  al.,  1994;  Lashlee,  1994).  Amongst  these 
robust  control  techniques,  mixed  norm  H2/H„o  optimal  control  methods  have  many 
advantages  for  designing  controllers  for  smart  structures.  These  control  techniques 
combine  the  good  closed  loop  performance  of  H2  controllers  with  the  good  stability 
robustness  of  controllers.  Mixed  norm  H2/H^  controllers  are  less  conservative  for 

'Department  of  Electrical  Engineering  and  Intelligent  Systems  Center.  320  ERL  Bldg,,  University  of  Missouri-Rolla, 
Rolla,  MO  65401,  U.S.A. 

'General  Motors  Corporation,  Warren,  Ml,  U.S.A. 
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systems  with  known  disturbance  power  spectral  densities  and  poorly  modeled 
disturbances  and  provide  a  trade-off  between  H2  and  performance.  Three  different 
approaches  for  the  formulation  and  solutions  of  H2/H^  controllers  are  available  in  the 
literature  (Bernstein  and  Haddad,  1989;  Bernstein,  Haddad  and  Nett,  1989;  Doyle,  Zhou 
and  Bodenheimer,  1989;  Khargonekar  and  Rotea,  1991).  The  application  of  H2/H„ 
controller  design  methodologies  for  smart  structures  is  reported  in  this  paper.  We  also 
propose  to  develop  distributed  sensors  for  the  measurement  of  physical  quantities  of 
structural  systems.  These  measured  quantities  can  be  used  for  the  implementation  of  H-,/ 
Hoo  controllers. 

The  smart  material  polyvinylidene  fluoride  (PVDF)  film  has  been  utilized  in  the 
development  of  customized  distributed  sensors  for  simple  cantilever  beams.  The.se 
measurements  are  used  for  the  identification  of  mathematical  models  of  the  structure 
from  experimental  data.  A  simple  direct  state  feedback  controller  has  been  demonstrated 
using  the  proposed  distributed  sensors. 


MIXED  H2/H^  optimal  CONTROL  TECHNIQUES 


Recently  there  has  been  a  great  deal  of  interest  in  the  formulation  and  solution  of  a 
mixed  H2/H^  control  methodology  which  can  handle  bounded  spectrum  and  bounded 
energy  inputs  simultaneously.  These  methods  promise  robustness  and  performance  in  the 
controller  design. 


CO, 

0)2 


Zl 

Z2 

y 


Figure  1.  Standard  Mixed  Norm  H2/H^  Problem 


The  design  objective  is  to  find  an  admissible  controller  K(s)  in  Figure  (1)  such  that 
(i)  ||z,||2  is  minimized  in  the  face  of  O),  while  at  the  same  time  (ii)  <  r  holds,  where 
Tj  is  the  transfer  function  between  0)2  to  Z2  -  The  infinity  norm  ||7'2||^<r  guarantees  a 
given  level  of  robustness,  and  the  smaller  value  or  r  yields  a  better  robust  system.  A 
trade  off  between  the  two  design  objectives  has  to  be  made  in  order  to  satisfy  both  of 
these. 


A.  BERNSTEIN  AND  HADDAD’S  ALGORITHM 

Among  the  three  groups,  Bernstein  and  Haddad  were  the  earliest  to  work  on  the 
mixed  norm  H2/He^  problem  (Bernstein  and  Haddad,  1989;  Bernststein,  1989).  In  their 
work,  a  special  case  of  the  general  mixed  norm  H2/H^  problem  was  considered,  where 
only  one  exogenous  input  signal  was  used.  That  is,  co,  =  0)2  =  to,  where  the  disturbance 
set  to  is  interpreted  as  a  standard  white  noise  signal. 

Correspondingly,  T^  and  are  defined  as  transfer  functions  from  O)  to  z,  and  to  to 
Z2  respectively.  The  nth  order  plant  P  in  state  space  form  is  given  by 
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A  D,  B 

p  —  El  0  E2 
•El  00  0  ^2^ 

C  Dj  0 

Alternatively,  the  problem  can  be  formulated  a  follows; 

Given  a  nth  order  plant 

X  =  Ax(t)  +Bu(t)  +D|CO(0 
y  =  Cxit)  +020(1) 

a  nth  order  dynamic  compensator  given  by 

u(t)  =  C,x,(t) 

can  be  determined  such  that  it  satisfies  the  performance  criterion 


(1) 


(2) 


(3) 


J(A^,B^,CJ  =  limE  {x^RiX  +  u^R2u}  (4) 

f  — >  00 

with  an  H^o  performance  bound  on  the  closed  loop  transfer  function  H  (s)  from  ©  (t)  to 
Ei^x(t)  +E2„uit)  given  by 


I|//(^)IL^Y  (5) 

The  compensator  satisfying  the  conditions  in  Equations  (4)  and  (5)  is  given  by 

A,  =  A-Qi:-'LPS  +  y^QRi„ 

=  QC'Vl'  (6) 

C,  =  -RI'B^PS 

where  Q ,  P  and  Q  satisfy  the  coupled  Riccati  equations 

0  =  AQ  +  QA^ +Vi+y^QRiM-Q^Q  (2) 

0=  (A  +  Y''[e  +  G]E,J  V  +  P(A  +  y‘'[(2  +  e]/?,J  +E,-5^PZP5  (8) 

0  =  (A  -  yP5  +  y’^GE,  J  Q  +  Q(A-  I.PS  +  y^QRiJ 
+  y"Q  (E,.  +  p'5VlP5)  Q  +  QLQ 


(9) 
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The  algorithm  described  above  assumes  the  original  system  is  a  strictly  proper 
system.  Usually  structural  systems  are  described  by  proper  transfer  functions  rather  than 
strictly  proper  systems.  Hence  the  Bernstein  and  Haddad  algorithm  has  been  modified  to 
accommodate  the  proper  transfer  function  system. 

If  we  remove  the  restriction  on  the  strictly  proper  system,  P  becomes 

A  Z),  B 
p  ^  E,  0  E, 

E\^  0 
C  D 

Using  the  following  solutions  of  the  Riccati  equations  for  Q,  Q  and  P : 

AQ  +  QA^  +  -QCL-i-^R,J)Q+V,  =  0 
(A+Y'"[(2  +  fi]/?,J^P  +  B(A  +  Y"'[e  +  e]B,J  +R,-S'PY.PS  =  0 


(10) 

(11) 

(12) 


0  =  (A  +  f^QR,^  -  ZP5)  e  +  e  (A  +  -  ^PS)  ^ 


+  Q'^Q  +  i^QRvM  + 

the  controller  (A^,  CJ  is  obtained  as 

A,  =  A-i:PS-Q'L  +  f^QR,^  +  QdVjDRl^B'PS  (14) 

Be  =  QC'r,'  (15) 

C,  =  -R-^B^PS  (16) 

and  the  cost  function  J  is  computed  by  substituting  q  and  R . 

J  =  tr{{Q  +  Q)R,  +  QS'pi.PS}  .  (17) 


B.  KHARGONEKAR  AND  ROTEA’S  ALGORITHM 

The  real  mixed  control  problem  (MCP)  was  considered  first  by  Rotea  and 
Khargonekar  in  (Khargonekar  and  Rotea,  1991)  towards  solving  a  special  case  of  full 
state  availability.  That  is,  no  upper  bound  to  the  two-norm  is  used  and  the  real  two-norm 
is  minimized.  Two  sets  of  exogenous  input  signals  and  two  sets  of  controlled  outputs  are 
allowed. 

Let  the  state-space  description  of  P  be  given  by 


X  =  Ax{t)  -I-  B|(0|  +82(^2  +  Bu 


(18) 
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z,  =  C,x  +  D|M 

(19) 

Z2  =  C2X  +  D2U 

(20) 

y  =  X. 

(21) 

Two  different  problems  are  considered  in  this  technique: 

Problem  A.  Constrained  optimal  control  problem:  For  the  plant  P ,  find  and  admissible 
(internally  stabilizable)  controller  K  that  achieves 

K  admissible  and  |7'2||^<  y} 

where 

r,  =  ^  (22) 

0), 

^  (23) 

02 

Problem  B.  Find  an  admissible  controller  K  that  achieves 

/n/{(|r,||2:  K  admissible} 

and  such  that  ||7’2||„  <  y. 

Note  that  Problem  A  represents  a  constrained  optimization  problem  (with  the 
constraint  ||7’2||„<y),  whereas  Problem  B  represents  an  unconstrained  problem  of 
minimizing  an  H2'porformance  measure  that  also  satisfies  an  Hoo^norm  bound.  If  an 
admissible  K(s)  exists  and  solves  Problem  B,  then  this  K{s)  is  also  a  solution  to 
Problem  A.  Thus  Problem  B  provides  sufficient  conditions  for  the  existence  of  a  solution 
to  Problem  A.  The  solution  method  minimizes  an  upper  bound  of  the  two-norm,  but  not 
the  two-norm  itself.  The  solution  of  this  problem  requires  full  state  feedback.  The 
distributed  sensor  proposed  in  the  next  section  will  provide  the  measurement  of  states  for 
the  implementation  of  H2/H^  control. 

DISTRIBUTED  SENSORS  FOR  STATE  MEASUREMENT 

The  measurement  of  the  model  states  of  a  structural  systems  is  possible  by  utilizing 
appropriately  shaped  PVDF  film.  PVDF  film  is  a  piezoelectric  material  available  in  thin 
sheets  which  can  easily  be  cut  to  the  desired  shape.  The  charge  developed  on  PVDF  film 
is  given  by 

qit)  =  Kf{b{x) — -y(x,t)dx  (24) 

,,  3X 

where  A}  is  a  constant  which  is  a  function  of  the  film  and  structure  parameters  and  b  (x) 
is  the  shape  of  the  sensor.  The  current  generated  by  the  film  is  proportional  to  the 
derivative  of  the  charge  since  the  film  is  a  capacitive  device.  For  the  cantilever  beam 
investigated  in  this  paper,  the  parameters  we  are  interested  in  measuring  are  the 
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displacement,  rotation,  velocity  and  angular  velocity  and  discrete  points  on  the  structure. 
The  shapes  required  for  the  point  measurement  of  these  quantities  are  derived  in  (Butler 
and  Rao,  1994)  and  are  very  simple  geometric  shapes.  This  allows  for  the 
implementation  of  several  sensors  with  a  single  sheet  of  film. 

The  measurement  of  point  displacement  is  accomplished  by  measuring  the  charge 
from  a  triangularly  shaped  PVDF  film  sensor.  The  sensor  segment  labeled  S2  in  Figure 
(2)  demonstrates  the  shape  of  the  sensor  for  the  measurement  of  the  displacement  of  the 
beam  at  position  L/2.  The  geometric  sum  of  sensor  segments  S2,  S3  and  S5  show  the 
shape  for  the  measurement  of  the  tip  displacement.  The  rotation  at  a  point  on  the 
cantilever  beam  is  measured  by  taking  the  charge  from  a  rectangularly  shaped  PVDF 
film.  Sensor  segments  SI,  S2  and  S3  show  the  shape  of  the  sensor  for  the  measurement 
of  the  rotation  of  the  beam  at  L/2.  The  summation  of  the  charge  from  all  five  of  the 
segments  results  in  the  rotation  of  the  tip  of  the  structure. 


Figure  2.  Segmented  PVDF  Film  on  Cantilever  Beam 


The  measurement  of  the  point  velocity  of  the  structure  is  accomplished  by 
measuring  the  current  on  a  PVDF  film  sensor  with  the  same  shape  as  the  displacement 
sensor  described  above.  The  measurement  of  the  point  rotation  rate  is  made  be 
measuring  the  current  from  a  sensor  with  the  same  shape  as  the  rotation  sensor  described 
above.  As  shown  in  Figure  (2),  several  sensor  shape  functions  can  be  obtained  by 
partitioning  a  single  sheet  of  PVDF  film  into  electrically  isolated  segments.  A  summary 
of  the  quantities  that  can  be  measured  with  the  segmented  PVDF  sensor  shown  in  Figure 
(2)  are  listed  in  Table  I. 


TABLE  1.  Sensor  Signals  in  Figure  (2) 


Physical  Parameter 

Segment  Numbers 

Measurement  T^pe 

Displacement  at  L/2 

S2 

Charge 

Rotation  at  L/2 

S1+S2+S3 

Charge 

Velocity  at  L/2 

S2 

Current 

Angular  Velocity  at  L/2 

S1+S2+S3 

Current 

Displacement  at  L 

S2+S3+S5 

Charge 

Rotation  at  L 

S1+S2+S3+S4+S5 

Charge 

Velocity  at  L 

S2+S3+S5 

Current 

Angular  Velocity  at  L 

S1+S2+S3+S4+S5 

Current 
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SYSTEM  IDENTIFICATION  USING  POINT  MEASUREMENTS 

In  addition  to  simplifying  controller  implementation,  the  point  measurements 
described  in  the  previous  section  aid  in  the  experimental  system  identification  of  the 
structure.  This  section  describes  a  method  which  has  been  developed  for  the 
experimental  identification  of  a  state  space  model  of  a  cantilever  beam  given  by 

‘  X  =  Ax  +  Bu 

.  (25) 

y  =  Cx 

Models  resulting  from  this  system  identification  procedure  closely  resemble  the  finite 
element  model  of  the  structure  and  have  states  with  the  same  physical  significance.  The 
states  of  the  experimental  model  are  displacements,  rotations,  velocities  and  angular 
velocities.  Unlike  the  finite  element  model,  it  is  not  necessaries  to  generate  a  large  model 
(several  elements)  in  order  to  accurately  represent  a  few  modes. 

Identification  of  the  cantilever  beam  model  begins  by  measuring  the  eigenvalues  of 
the  system  that  we  are  interested  in  modeling.  In  general,  this  can  be  done  with  swept 
sine  measurements  or  by  simply  tuning  the  sinusoidal  input  frequency  of  the  actuation 
signal  to  achieve  the  maximum  observed  response  on  any  of  the  sensors.  With  the 
frequencies  co^  known,  the  damping  of  the  mode  can  be  measured  by  turning  off  the 
actuator  once  the  steady  state  response  is  achieved  and  measuring  the  decay  of  the  time 
domain  signal.  Performing  these  measurements  for  each  of  the  modes  of  interest  allows 
for  the  calculation  of  the  matrix  A  given  by 

A  =  r'AT,B  =  r'B,C  =  CT  (26) 

where  A  is  the  diagonalized  system  matrix.  The  linear  transformation  in  Equation  (26) 
given  by  r  =  [p,  P„  consists  of  columns  of  eigenvectors  P,  of  the  system. 

The  eigenvectors  of  the  system  matrix  can  be  measured  using  a  sinusoidal  input 
signal  at  the  frequency  of  the  eigenvalue  and  measuring  the  steady  state  magnitude  and 
phase  of  each  of  the  distributed  sensors.  These  quantities  represent  the  eigenvectors  P, 
of  the  system  and  completely  define  the  linear  transformation  to  the  diagonal  system.  For 
the  cantilever  beam  with  no  significant  actuator  dynamics,  all  of  the  measurements  will 
be  in  phase  or  90°  out  of  phase  with  the  input  signal.  The  system  A  matrix  can  be 
calculated  by 

A  =  tat'  (27) 

The  B  matrix  in  Equation  (26)  of  the  diagonal  system  has  the  form 


B  =  [bj 


-bj  bJ 


so  that  the  state  equation  for  a  fourth  order  diagonal  system  is 

[bJ  1 


X,  =  Ax,  +  BUj  =  Ax,  +  w, 

bJ 
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where  m,  is  the  input  corresponding  to  the  ith  mode.  The  elements  of  the  diagonal 
matrix  B  can  be  calculated  from  the  data  collected  for  the  identification  of  the  A  matrix. 
Solving  Equation  (29)  for  B  with  the  steady  state  sinusoidal  input  m,  =  q,cos  (o)/) 
gives  ^ 

,  X:  -  AXi  7^' JC,  -  T'AX: 

bi  =  - - ■'  =  - ^ (30) 


Since  the  vector  x  consists  of  sinusoids,  x  can  be  calculated.  The  time  derivative  of  the 
measured  states  is  given  by 


X,.  =  (oxj 


yielding 

<li 

Solving  Equation  (32)  for  each  of  the  individual  eigenvalue  inputs  gives  the  B  matrix. 
The  B  matrix  can  be  determined  using  Equation  (26). 

Since  the  physical  significance  of  the  states  is  well  known  for  this  system,  the  C 
matrix  can  be  selected  as  desired. 


CANTILEVER  BEAM  FULL  STATE  FEEDBACK  DESIGN  EXAMPLE 

The  measurement  of  the  states  of  a  structural  system  allows  for  the  implementation 
of  full  state  feedback  controllers  such  as  pole  placement,  LQR,  Ho,,  and  H2/H„„  controllers 
using  only  simple  analog  gain  circuits.  The  linear  quadratic  regulator  (LQR)  controller 
has  good  robustness  properties  and  will  be  used  here  to  illustrate  the  simplicity  of  the 
full  state  controller  implementation.  The  measurement  of  the  states  of  the  system 
eliminates  the  need  for  state  estimation.  This  means  that  the  robustness  properties 
inherent  in  the  LQR  controller  are  retained  in  the  controller  implementation  and  do  not 
need  to  be  recovered  with  a  loop  transfer  recovery  methodology. 

A  finite  element  model  of  the  structure  shown  in  Figure  (3)  will  be  utilized  to 
illustrate  the  benefits  of  the  controller.  As  mentioned  in  the  previous  section,  an 
experimental  model  can  be  determined  for  the  actual  system  which  accurately  models  all 
of  the  modes  present  in  the  model.  The  error  inherent  in  finite  element  model  using  only 
one  segment  would  not  be  present  in  the  experimentally  determined  model. 


Figure  3.  Cantilever  Beam  with  Segmented  PVDF  Film  Sensor 
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The  distributed  sensor  shown  in  Figure  (3)  will  provide  the  measurement  of  the  tip 
displacement  (charge  on  SI),  rotation  (charge  on  S1+S2),  velocity  (current  from  SI)  and 
angular  velocity  (current  from  S1+S2).  The  mass  and  stiffness  matrix  for  a 
corresponding  single  element  finite  element  model  of  the  beam  is  given  by 


1.132e-2  -4.459e-4 

3.261e2 

-4.556el 

(33) 

__ 4.459^-4  2.265e-5j’ 

_-4.556el 

8.486e0_ 

The  FEM  state  space  model  with  1%  damping  on  each  of  the  modes  is  given  by 

l.Olel  -4.18el  leO  0  0  1  0  0  0  0 

5.62el  -1.24el  0  leO  0  (34) 

2.25e5  -4.78e4  l.OleO  -4.18e-l  ’  -5.18e2  ’  0  0  1  0  ’  0 

6.44e6 -1.32e6  5.62el  -1.25elJ  [-1.71mJ  [o  0  0  ij  [Q. 

In  the  LQR  controller,  the  cost  function 


J  =  j  Qx  +  Ru)  dt  (35) 

0 

is  minimized  subject  to  the  constraint  equation  x  =  Ax  +  Bu  The  cost  on  the  states  and 
the  cost  on  the  control  signal,  Q  and  R  respectively,  were  chosen  to  provide  11% 
damping  on  each  of  the  modes. 


2.11  el  -4.03e6  0  0 

-4.03e6  8.08e5  0  0 

0  0  1.50e2  -8.51e0 

0  0  -8.51e0  9.23e-l 


R  =  30000 


resulting  in  the  feedback  gain  matrix 

K  =  [9.3 e-1  -1.73e-l  1.05e-l  -8.57 e-3] 
Simulations  of  the  open  and  closed  loop  response  is  shown  in  Figure  (4). 
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Figure  4.  Simulated  Open  and  Closed  Loop  Response  of  Cantilever  Beam 
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CONCLUSIONS 

A  method  to  utilize  multiobjective  robust  controllers  for  smart  structures  is 
presented  in  this  paper.  The  proposed  method  provides  a  trade-off  between  closed-loop 
performance  and  disturbance  attenuation.  A  simple  modification  is  provided  for 
designing  controllers  for  systems  with  proper  transfer  functions. 

We  have  also  developed  distributed  sensor  shape  functions  for  the  measurement  of 
physical  quantities  of  structural  systems  using  PVDF  film  sensors.  These  sensor 
measurements  can  be  used  for  the  implementation  of  full  state  feedback  H2/H00 
controllers  on  smart  structures.  We  have  demonstrated  the  applicability  of  distributed 
sensors  for  the  implementation  of  full  state  feedback  systems.  We  are  currently  working 
on  implementation  of  these  algorithms  on  test  articles.  We  have  also  successfully 
demonstrated  a  procedure  for  system  identification  using  point  measurements. 
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